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Abstract 

We study the evolution of the interface given by two incompressible fluids with different 
densities in the porous strip R x [—1, 1}. This problem is known as the Muskat problem and 
f*"*- ■ is analogous to the two phase Hele-Shaw cell. The main goal of this paper is to compare the 

qualitative properties between the model when the fluids move without boundaries and the 
model when the fluids are confined. We find that, in a precise sense, the boundaries decrease 
| the diffusion rate and the system becomes more singular. 
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1 Introduction 

In this paper we study the evolution of the interface between two different incompressible fluids 
with the same viscosity in a flat two-dimensional strip. This problem has an interest because it 
is a model of an aquifer or an oil well, see [18]. In this phenomena, the velocity of a fluid in a 
porous medium satisfies Darcy's law 

rS [ 

C3; -v = -Vp-gpe 2 , (1) 

K 

where p is the dynamic viscosity, k is the permeability of the medium, g is the acceleration due 
to gravity, p is the density of the fluid, p is the pressure of the fluid and v is the incompressible 
field of velocities, see [3, 19]. 

Equation (1) has also been considered as a model of the velocity for cells in tumor growth, 
see for instance [13, 221 and references therein. 
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Figure 1: Physical situation for an interface z(a,t) in the strip R x (—1,1). 



The motion of a fluid in a two-dimensional porous medium is analogous to the Hele-Shaw 
cell problem, see [14]. In this case the fluid is trapped between two parallel plates. The mean 
velocity in the cell is described by 

— v = -Vp - gpe 2 , 

where b is the (small) distance between the plates. 

We consider the two-dimensional flat strip 5 = lx(-I,l)cK 2 with I > 0. In this strip we 
have two immiscible and incompressible fluids with the same viscosity and different densities, p 1 
in S' 1 (t) and p 2 in S 2 (t), where S l (t) denotes the domain occupied by the i— th fluid. The curve 

z(a,t) = {(zi(a,t), Z2(a,t)) : a G R} 

is the interface between the fluids. We suppose that the initial interface fo(x) is a graph and 
l/ol^)! < I for an x - The character of being a graph is preserved at least for a short time (see 
Section 3). The Rayleigh- Taylor condition is defined as 

RT(a,t) = -(Vp 2 (z(a,t)) - Vp l (z(a,t))) ■ d^z(a,t). 

Due to the incompressiblity of the fluids and using that the curve can be parametrized as a 
graph, the Rayleigh- Taylor condition reduces to the sign of the jump in the density: 

RT = g{p 2 -p 1 ) >0. 

This condition is satisfied if the denser fluid is below. 

We consider the velocity field v, the pressure p and the density p 

P(t) = p 1 XsHt) + p 2 XsHt), (2) 

in the whole domain S. We also consider the conservation of mass equation, so we have a weak 
solution to the following system of equations 

-v(x, y, t) = -Vp(x, y, t) - gpe 2 in S, t > 0, 

< V-v(x,y,t) = in S, t > 0, (3) 

dtp(x, y,t) + v ■ Vp(x, y,t) = in S, t > 0, 
k /(x,0) = / o (a:) inR, 
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with impermeable boundary conditions for the velocity (see Section 2). 

We denote by v 1 (x, y, t) the velocity field in and by v 2 (x, y, t) the velocity field in S 2 (t). 
Because of the incompressibility condition, the normal components of the velocities f 1 ,^ 2 are 
continuous through the interface. Moreover, the interface moves along with the fluids. Therefore, 
if initially we have an interface which is the graph of a function, we have the following equation 
for the interface: 



d t f(x,t) = (-d x f(x,t),l)-v i (x,f(x,t),t) = y/l + {d x f{x,t)fn-v\ (4) 

where n denotes the unit normal to the interface. 

In each subdomain S l {t) the fluids satisfy Darcy's law (1), 

^v i (x,y,t) = -Vp l (x,y,t) -gp%l) in S l (t), (5) 

K 

and the incompressibility condition 

V-v i (x,y,t) = Q in 5*'(t). (6) 
We define the following dimensionless parameter (see [4] and references therein) 

This parameter is called the nonlinearity (or amplitude) parameter and we have < A < 1. 

The case A = 1 is the case where / reaches the boundaries and we call it the large amplitude 
regime. In [ll ] they consider a two dimensional droplet in vacuum over a plate driven by surface 
tension. 

The case A = is the deep water regime for which the equation reduces to 

BJ = ^P.V. / ^[fj-^f-^ d,,. (8) 
2vr J R rj 2 + (/(s) - f(x - r?)) 2 

It has been shown, for equation (8), local existence in Sobolev spaces when the Rayleigh- Taylor 
condition holds (see [9]), a maximum principle for the L°° norm of / and also a maximum 
principle for (see [10]). For initial data with [I^/oIIl 00 < 1 follows global existence 

of W 1 ' 00 solution (see [i]). For large initial datum there are turning waves, i.e a blow up for 
(see [6]). For other results see [1, 5, 6, 8, 15, 23]. 
The equation for the evolution of the interface in our bounded domain, which is deduced in 
Section 2, is 

2 _ „1 

(d x f (x) - d x f (x - rj)) Si (x, 77, /) 



+ {d x f (x) + d x f (x - 7]) E 2 (x, rj, f) 



dv = ^^A[f)(x), (9) 



where the singular kernels Si and S2 are defined as 



" sinh (f 77) 

V) = cosh (§r,) - cos(§(/(x) - f(x - v))) ' (10) 



corresponding to the singular character of the problem, and 

sinh (§77) 



V) cosh (f lV ) + cos(§(/(x) + /(* - 77))) ' (U) 
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which becomes singular when / reaches the boundaries. The text P.V. denotes principal value. 
As for the whole plane case (see [7, 11]) the spatial operator can be written as an 

x-derivative. Indeed, 



tan( ^/(EKfa) 



A\f](x) = -P.V. / d x arctan ^ — - 

1 n ' * A ^ ^ tanh(§2) 

+ § P.V. ^ (-an ^ (|f))) 



di] 

■ — nl \ / 7r n \ \ \ 

dr? (12) 



and we conclude the mean conservation 



/ f(x,t)dx= [ f (x)dx. (13) 
Jr Jr 

When we do not parametrize the curve as a graph, i.e., we consider z(a) = (21(a), 22(a)), 
we obtain the equation 



d t z = ^P.V. 



+ 



(d a z(a) - d a z(r])) smh(zi(a) - 21(77)) _ 

) 

dr). (14) 



cosh(zi(a) - 2:1(77)) - cos(z 2 (a) - ^(v)) 
(d a z 1 (a) - d a z 1 (7]),d a z 2 (a) + d a z 2 (rj)) sinh(zi(a) - 2:1(77)) 



cosh(zi(a) - zi(rj)) + cos(z 2 (a) + z 2 (r])) 
As our interface moves in a bounded medium the correct space to consider is 

Hf = H s (R)n{f: H/llxoo </}. 

The density p defined as in (2) is a weak solution of the conservation of mass equation present 
in (3) if and only if the interface verifies the equation (4) (see Proposition 1 in Section 2 below). 
It also follows (see 21) that if we take the limit A — > we recover the equation (8) (see [9]). 

In a recent work [12], J. Escher and B-V.Matioc studied the problem (5), (6) in the case 
with different viscosities and surface tension in a periodic (in x) domain when < A < 1. 
They obtained an abstract evolution equation for the interface and showed well-posedness in 
the classical sense when the Rayleigh- Taylor condition is satisfied and the interface is in a 
neighbourhood of the zero function in certain Holder spaces. They consider the problem as a 
problem in two coupled domains. The domains are coupled by the interface and by the Laplace- 
Young condition 



p 2 \ 



x, f{x, t),t) - p 1 (x, f(x, t),t) = 7k[/], 



where n[f] denotes the curvature of the interface f(x,t) and 7 denotes the surface tension 
coefficient. 

In Section 3 we study the similarities between the case .4 = and < A < 1; first, 
we prove local well-posedness in Sobolev spaces (see Section 3.1) and instant analyticity in a 
growing complex strip when the Rayleigh- Taylor condition is achieved (see Section 3.2). The 
last similarity studied in Section 3.3 is that for arbitrary initial curves which are analytic there 
is an unique local solution, which is analytic, both forward and backward in time. We remark 
that for this result the Rayleigh- Taylor condition is not needed. The proofs follows the steps 
of the paper [6] . Here we show that the contribution from the boundary does not affect the a 
priori estimates from [6]. 

The main purpose of this work is to study the differences between the case with infinite depth 
and the case with bounded medium. This is done in Section 4, where we study some qualitative 
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properties of the solution. We prove the maximum principle for ||/(t)||z<» and also for ||c?a;/ ||z,°° 
by studying the evolution of the maximum or the minimum values. The ODEs coming from 
these analysis have local and non-local terms and the main dificulty is to compare these two 
different kind of terms in order to ensure the decay. In particular, we prove the following decay 
estimate for ||/(£)||z,°o in a confined medium: 

^H/(*)IU~ <-c(\\fo\\^A\foh^,P^p\l)e~™^. (15) 

As a corollary we conclude that the unique one-signed, integrable, stationary solution is the rest 
state. Let us observe that the natural boundary condition for the velocity, v ■ n = 0, imposes 
that if our initial interface is close enough to the boundary the evolution of the maximum is 
very slow. Due to this fact, we obtain the slow decay inequality (15). 

We show that if the initial data is in a region depending on ||/o||l°° and ||(? x /o||l°° then we 
have a maximum principle for ||c?x/||_L°° i n a confined medium. We consider a smooth initial 
data /o in the Rayleigh- Taylor stable regime such that the following conditions holds: 

\\d t foh« < 1, (16) 
t ail < \\d x f \\Loo tanh (£) , (17) 



and 



\d x fo\\L~ + |2(cos [-) - 2)sec 4 [-) \\\d x f \\U) ^ 



(u taa(^ l|a " / °"^ 00 



tanh(f) 




6 tanh (§) 4Z 2 



+ 4tan (^H/oIIl-) - ^l^/olk- cos (j||/o|U«) < (18) 



Then 



\\d x f{t)\\»° < \\dxh\\L°°- (19) 
Moreover, if (x(l),y(l)) is the solution of the system 

tan(ff) -j/tanh(§) =0 



i+!,[!l+t ^M)) ( 2U, 



(y + |2(cos (§) - 2) sec 4 (f ) |y») V ^£gf^ (§) 

+4 tan (§x) - 4ycos (fx) = 0, 

and we have that HSx/olli 00 < 2/(0 an d H/olli 00 < x(l) we have that 

R/Hloo < 1. 

The effect of the boundaries is very important at this level, and we obtain a region (the 
yellow one in Figure 2) where we do not have maximum principle for ||3 x /o|| but we have an 
uniform bound 11^/(4)11^00 < 1 Vt > 0. The green region is the region where we have maximum 
principle for the derivative. Due to the term coming from the effect of the boundaries, H2, 
the conditions that we obtain are much more restrictives than in the deep water regime (the 
case with infinite depth). The previous result gives us conditions on the smallness of A and 
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Figure 2: Different regions in (||/o||l°°, HSe/oIIz, 00 ) for the behaviour of when tt = 21. 

II^x/oIIl 00 (which, for fixed amplitude, can be understood as the 'wavelength' of the wave) so, 
roughly speaking the Theorem says that if we are in the long wave regime (small amplitude and 
large wavelenght) then there is no turning effect, i.e. there are no shocks. We remark that if we 
take the limit A — > we recover the result for the deep water regime contained in [10] 

We study the formation of singularities in Section 4.2. The singularity is a blow up of 
||9e/(*)||l°°- Physically this result means that there are waves such that they 'turn over'. 
Moreover, we compare this result with the result for the deep water regime (see [(>]). In particular, 
we prove the following turning effect in a confined medium: There exist initial data Zq{ol) = 
(z\(a), Z2(a)) such that in finite time the solution of (14) achieves the unstable regime only 
when the depth is finite. If the depth is infinite the same curves are depleted. 

In Section 4.3 we show an ill-posedness result in Sobolev spaces. The key point of this result 
is that we do not need global existence for some class of solutions to prove the result (compare 
with [ ] and [23]). 

Remark 1 In order to simplify the notation we take p,/n = g = 1 and we sometimes suppress 
the dependence on t. We denote Vi the component i— th of the vector v. We remark that v % is 
the velocity field in S l {t). We write n for the unitary normal to the curve T vector and n for 

2 1 

the non-unitary normal vector. We denote p — — ^ — . 

2 The equation for the internal wave 

In this section we obtain the equation for the interface z(a, t) in an explicit formula. First we 
have to add impermeable boundary conditions for v, i.e. v(x, ±Z, t) ■ n = 0. 

Using the incompressibility condition we have that there exists a scalar function ^ such that 
v = V J ~ 1 I / . The function ^ is the stream function. Then 



= -curl(0, p) = u) 



where the vorticity is supported on the curve 



u(x,y) 



w(a)S((x,y) - z(a,t)) 



with amplitude 



ro(a) 



(p 2 - p 1 )d a z 2 (a). 
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In this domain we need to obtain the Biot-Savart law. The Green function for the equation Au = 
f in the strip K x (0, 21) (with homogeneous Dirichlet conditions) is given by the convolution 
with the kernel 

1 

G{x,y,n,v) = — log ( y/{x - fi) 2 + (y- (4nZ + z/)) 2 ) 

- log (y/{x - /") 2 + {y~ (4n/ - v)) 2 ) 



n=— oo ■- 



The Biot-Savart law in this strip is given by the kernel 



BS{x,y,n,v) = V^ y G{x,y,n,v) = — ^ 



(r+) x (r-)- 
K\ 2 



I — 12 

J- ra 



where 



r+ = (x-/i,j/-(4nZ + i/)), T n =(x-n,y-(4nl-v)). 
It is useful to consider complex variables notation. Then 



^ oo 

BS(x,y,n,v) = — Y 



n=— oo <- 



1 1 

1 n L n 



Fixed n, we compute the following 

1 1 



2{x - fj, + i(y - v)) 



r+ r± n (x- / u + %-^)) 2 + (4n0 2 ' 



1 1 

F + r~ 

1 n 1 — j 



2(x - fi + i(y + u)) 



(x — fjb + i(y + v)) 2 + (4nl) 2 ' 

We change variables (y — I = y, v — I = v) to recover the initial strip S 
without lossing generality we take I = n/2. Due to the formula 



x (— I, I), moreover, 



2z 



— coth ( - ) , 



oo 
n=l 

we obtain that the Biot-Savart law in cartesian coordinates is given by 



BS{x,y,n,v) = — I 



sin(y — z/) 



sin(y + v) 



Air \ cosh(x — fi) — cos(y — v) cosh(x — fi) + cos(y + v) ' 

sinh(x — fi) sinh(x — /j,) 



cosh(x — n) — cos(y — v) cosh(x — fi) + cos(y + v) J ' 
Using the formula for the vorticity we have that the velocity is 



v(x,y) = — j -(,]) 



sin(y - z 2 (f3)) 



+ 



sin(y + z 2 (/3)) 



cosh(x - zi((3)) - cos(y - z 2 (f3)) cosh(x - z\((3)) + cos(y + z 2 (f3)) ' 
sinh(x — z±(f3)) sinh(x — z\{j3)) 



cosh(x - zi{(3)) - cos(y - z 2 (fi)) cosh(x - z\((3)) + cos(y + z 2 (f3)) 
We use the identity 

<9,,log(cosh(zi(a) - 21(77)) ±cos(z 2 (a) ± z 2 (rf))) = 



d/3. 
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to obtain that the average velocity in the curve is 



v{z(a)) 



P-, 



P.V. / d aZl (ri) 



sinh(zi(a) — 21(77)) 



+ 



cosh(zi(a) — z\(r])) — cos(z2{a) 
sinh(zi(a) — 21(77)) 



cosh(zi(a) 
da ^2(77) 



21(77)) + cos(z 2 (a) + 22(77)) 
sinh(zi(a) — 21(77)) 



^2(7?)) 
d-q, 



cosh(zi(a) - 21(77)) - cos(z 2 (a) - 22(77)) 
sinh(zi(a) — 21(77)) 



drj 



cosh(zi(a) - 21(77)) + cos (22 (a) + 22(77)) 

The interface is convected by this velocity but we can add any velocity in the tangential 
direction without altering the shape of the curve. The tangential velocity in a curve only changes 
the parametrization. We consider then the following equation with the redefined velocity 

dtz(a) = v(z(a)) + c(a)d a z(a), 

where 

sinh(2i(a) — 21(77)) 



cm 



£p.v. 

2 



cosh(2i(a) — 21(77)) 



cos (22 (a) 
+ 



22(77)) 



sinh(2i(a) — 21(77)) 



-drj. 



cosh(2i(a) - 21(77)) + cos(2 2 (a) + 22(77)) 

Following this approach we obtain (14). Because of that choice of c(a) we obtain that, if initially 
the curve can be parametrized as a graph, i.e., z(x,0) = (x,/o(x)), we have that the velocity 
v\ on the curve is zero, thus our curve is parametrized as a graph for t > and we recover the 
contour equation (9). 

Note that when I — > 00 in the equation (9) we recover the equation for the whole plane (8): 



lim jA[f](x) 



-P.V. 



7T 



v(9 x f(x) - d x f(x - 7/)) 



dr), 



(21) 



77 2 + (f(x) - f(x - 7])) 2 

where A[f] is the operator defined in (9). 

Furthermore, we obtain the pressure p (up to a constant) solving the equation 

-Ap = gdyp, 

with Neumann boundary conditions 

d n p\ y =i = -9P 1 , d n p\ y= _i = gp 2 . 

In this way we obtain v,p satisfying Darcy's Law and the incompressibility condition. It is 
easy to check that p(x, y, t) is a weak solution of the conservation of mass equation. 

Definition 1. Let v be an incompressible field of velocities following Darcy's Law. We define 
the weak solution of the conservation of mass equation present in (3) as a function satisfying 



T 



p(x, y, t)d t 4>(x, y, t) + v(x, y, t)p(x, y, t)V x , y (f>(x, y, t)dydxdt = 



for all (f) G C C °°(]R x (-1,1) x (0,T)). 

We conclude this section with the following result. 

Proposition 1. Let p be the function defined in (2). Then p is a weak solution of the conser- 
vation of mass equation (see Definition 1) if and only if f is a solution of (9). 

The proofs of these two results are straightforward and, for the sake of brevity, we left them 
for the interested reader. 
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3 Similar results between the two regimes 



In this section we show a group of results for < A < 1 similar to those in the regime A = 0. 
The proofs follow the same ideas but, due to the structure in our equation (9), with a second 
term coming from the boundaries present in our model, there are some difficulties. We show the 
well-posedness in Sobolev spaces when the Rayleigh- Taylor condition is satified, i.e. the denser 
fluid is below the lighter one. In the case where the Rayleigh- Taylor condition is not satisfied but 
our initial data is analytic we also have a well-posedness result by means of a Cauchy-Kovalevski 
Theorem. We also prove the smoothing effect of the spatial operator in (9), i.e. the solution 
becomes instantly analytic. 



3.1 Well-posedness in Sobolev spaces 

In this section we sketch the proof of local well-posedness in Sobolev spaces in the Rayleigh- 
Taylor stable case: 

Without loss of generality we take 21 = it and p = 2. We indicate the constants with a 
dependency on I as c(l). The proof follows the same lines as in [9]. In order to deal with 
the kernel Eg, the kernel corresponding to the effect of the boundaries, we define the following 
energy: 

E[m) = \\f\\ 2 H s(t) + \\d[f]\\ L o C (t), (22) 
where d[f] :R 2 xMMl + is defined as 

d[f](x,rj,t) = — — ] rr. (23) 

cosh(r/) + cos( j (x) + / (x — r])) 

The function (23) measures the distance between / and the top and floor ±Z (recall that, without 
losing generality, we are supposing that I = |). In other words, < oo implies that 

||/||l°° < f ■ So this is the natural 'energy' associated to the space Hf(W). We obtain 'a priori' 
energy estimates as in [9]: The integrals corresponding to the kernel Si are the more singular 
terms and can be bounded as in [9] because has a singularity with the same order. The integrals 
corresponding to the kernel H2 are harmless and can be bounded using the definition of d[f\. 
For instance, one of the integrals arising in the study of the third derivative, after an integration 
by parts, is 



\dlf{x) 



/ +/ ) d x E 2 (x,r))dr)dx = I in + I out , 



and we obtain 

sin (f(x) + f(x - r/)) (d x f(x) + d x f(x - rj)) 



|^/(x)| 2 P.V. / sinh(|77|) 

•75(0,1) 



(cosh (77) + cos((/(x) + f(x - rj)))Y 



drjdx 



< cmf\\cA\d 3 J\\h\\d[f]\\l^ < c(l)\\f\\U\d[f u2 



With these techniques we obtain 



2 f sinh(|r/|) ; n| ., 

±\\f\\ H3 < c (i m f] + lf. 
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In order to use classical energy methods we have to bound the evolution of \\d[f] in terms 
of the energy E[f]. With this method we need a bound on ||9t/||L°°- In order to do this we split 
dtf in two terms, one for each kernel 



d t f = A l + A 2 . 



We give the proof for the 



A l = P.V. / (d x f(x) - d x f(x - Ti^E^x, r,)d V . 



For the term corresponding to the second kernel, A2, the procedure is analogous. 
We split A 1 in its 'in' and 'out' parts, A 1 = A\ n + A° ut , with 



A?<c(i)\\d 2 j\\ L ~, 



and 



A out < 



P.V. / d x f(x)E 1 (x,r l )dr l 

S C (0,1) 



+ 



P.V. / -d x f(x - r])E 1 (x,r])dr] 

B C (0,1) 



We have that the integral 



P.V. 



/ 

JB C (0,1 



sinh(?7) 



(0,1) sinh ( j) 



Using this fact and the classical and hyperbolic trigonometric identities we can write 

sinh(?7) 



P.V. / a x /(a:)Hi(x, rfidq = d x f(x)P.V. [ 



(0,i) 2 sinh 2 (2) 



1 



1 



1 mnH{f(x)-f{x-r,))/2) \ + (d x f(x)) 2 
Xi sinh 2 (»;/2) 



drj. 



We define the function 



Q(v) 



sm 2 ((f(x)-f(x-r ] ))/2) 
sinh 2 (r?/2) 



We observe that Q(0) = lim T; ^o Qi 7 ]) = {9 x f(x)) 2 , and we can apply the Mean Value Theorem 
to this function Q in order to have the cancelation needed at infinity. We obtain 



PV. / d x f(x)E 1 (x,r])dri 

Bc(0,l) 



<c(l)\\d x f\\L°°(l + \\d x f\\ L ° 



To bound 



P.V./ -d x f(x - 7])Ei(x, rj)dr] = P.V. / d ri f(x-r])E 1 (x,7])dr] 
Jb c (o,i) Jb c (o,i) 

we integrate by parts. We conclude 

P.V./ -d x f(x -r))E 1 (x,i])dr) 
Jb c (o.i) 



<c(0H/||l<»(1 + ||S x /||loo). 



10 



Now we can prove the last estimate: We have that 

jd[f] = d[f} 2 sin(/(x) + f(x - r,)Mf(x) + d t f(x - rj)) < c(l)d[f] 2 \\d t f\\ 
Due to the previous bound for ||dt/||z,°° we obtain 

j t d[f]<c(l)d[f]\W]\\ME[f] + l) 2 , 

and integrating in time we conclude that 

d[f]{t + h) < d[f](t)e^ +h <md[f]\\ L oo(ElfKs)+irds_ 

Finally we have that 

d iUim .. \\d[f]\\L°°(t + h)-\\d\f]\\L~(t) r7Vprfl , n4 
at h->o h 

Putting all together we obtain the following bound 

^[/]<c(I)(£[/] + l) 5 , 

therefore 

E[f)(t) < - E[k] 



y-4E[f Q Mi)t + 1 

With this 'a priori' bound we can obtain the local existence and uniqueness of classical solutions 
using the energy method (see [17]). 

3.2 Smoothing effect 

In this section we sketch the proof of the instant analyticity for the classical solution (which 
exists due to the result in the previous Section) if we start with a smooth initial data. 

Without lossing generality we consider I = ir/2 and p = 2. The proof, as the one in [6], 
relies on some estimates for the complex extension of the function / on the boundary of the 
strip M = {x + i£, |£| < kt, V < t < T(/o)} for certain k, a constant that will be fixed later. 
We define 



2 



f \f{x + ikt)\ 2 dx+ f \f(x-ikt)\ 2 dx, \\f\\ 2 H3{B) = ||/||£ 2(B) + ||^/||| 2( 

cosh 2 (?y/3) 



d + [f](x + iC,v) 
<r\f](x + i£,ri) 



cosh(r/) + cos(/(x + i£) + f{x + i£ — rj)) ' 
sinh 2 (r//3) 

cosh(r/) — cos(/(x + i£) — f(x + i£ — rj)) ' 



and 



m(t)=minRe i + a 1 /(7))2 (24) 

We remark that as d[f] defined in (23), comes from the boundaries. d~[f] plays the role 

of the arc-chord condition (see [6]) and ensures that the singularity in the first kernel has order 
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two. In the case of complex functions we have that if the kernel Si is singular then the following 
equality holds 

cosh(r/) — cos(Re/(x + ikt) — Ref(x + ikt — 77)) cosh(Im/(x + ikt) — Imf(x + ikt — r/)) = 0. 
Assume now that |r?| > R > 21, where R >> 21 is a fixed constant, then 

- cosh(r/) + - cosh( J R) - cosh(2||/|| i0 o (B) ) < - cosh(r/) + - cosh(i?) - cosh(2||Im/|| i ^ (B) ) 
< cosh(r/) — cos(Re/(x + ikt) — Re/(x + ikt — rj)) cosh(Im/(x + ikt) — Im/(x + ikt — rj)), 



so we want that 



1 



cosh(i?) - cosh(2||/||_ L oo (B )) > 0, 



because in that case the kernel Hi is not singular. A similar analysis can be done for H 2 . We 
write 

Dim 



1 



cosh(i?) - 2cosh(2|/( 7 )|)' 

This term is a technical resource to obtain enough decay at infinity. We consider Hardy-Sobolev 
spaces (see [2] and references therein) on B so we want to obtain 'a priori' bounds on the 
following energy 



where 



Mf] = II/IIm) + \\d + [f]\\ L ~ m + ||<T[/]|| L ~(B) + P[/]|U~(B), 



\\F(x + i£,r])\\ L oo {M ) = sup \F(x + i£,rj)\. 



(25) 



The evolution for the complex extension of / is 

(d x f {x ± ikt) — d x f (x ± ikt — r/)) sinh (77) 



d t f(x ± ikt) = RV. 



cosh (if) — cos((/(x ± ikt) — f(x ± ikt — rj))) 

(d x f (x ± ikt) + d x f (x ± ikt - rj) sinh (rj) 
cosh (rj) + cos((f(x ± ikt) + f(x ± ikt — rj))) 

Recall that R >> 21 is a fixed constant. Then, as in [9, 6], we obtain 
^B[/]<exp( C (Z)(£ B [/] + l)) 



drj. (26) 



+ IIA 1 / 2 ^/!^ (k Im (V.V. jf Hi + H 2 ) 



+ 2k - 2vrm(t) , (27) 



where K = K(l) is an universal constant and /c is the width of the strip. At time t = 0, we have 
that if k is small enough (k only depends on the initial data) 



Ak+K 



Im P.V. 



-1 T =2 



t=0 



2^ 



L°°(B) 1 + II^e/oIIlo 



4fc- 



2tt 



l + \\d x fo\\L° 



< 0. (28) 



Then we need to show that this quantity remains negative (at least) for a short time. In order 
to do this we define the following new energy: 



<£b[/] = E v [f] + 



27rm(i) - K ||lm (P.V. / R Hi + H 2 ) | 



4k 



(29) 
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If 



< oo then < exp(c(l)(E^[f] + 1)) and we have the correct 'a priori' estimates. 

Now, if we have a classical 



We need to bound m'(t) and ^ ||lm (P.V. J R Hi + ^dr?) || ioo / 
solution with E®[f] < oo, the Sobolev embedding gives us that 



Re- 



G C^M x B) 



'l + (^/( 7 )) 2 

Thus we can apply Rademacher Theorem to m(t) in (24) and we get 

m / <exp(c(/)(€ B [/] + l)) 
Again, applying Rademacher Theorem, we have 



d_ 
~dt 



Im P.V. / Hi + E 2 d V 



<exp(c(0((£B [/] + !))■ 



We conclude 



and then we have 



-£b[/] < exp(c(Z) (<£ B [/] + !)), 



<£b [/](*) < -^log(exp(- C (0<S B [/o]) -c(0exp(c(0)i). 



Furthermore, there exists a time T = T(fo), T < T* = 
Now, for e > 0, we consider 



cx P (-c(0(i+e B [/o])) 

c(0 



e \ e . 

where J is the heat kernel, and the regularized problem 



dtp 



F e (t' 5 ), f e ' 5 (x,Q)=J e *fo(x) 



(30) 

such that <£ B [/] < c(/ )- 
(31) 

(32) 



where 



pes (je 



Je * P.V. 



+ J e *[ PV. 



(>7 e * d x f £ > s (x) - Je * d x f e > 5 (x - r))) sinh(7?) d 
cosh(?7) - cos(J e * / e ' 5 (x) - J € * / e ' a (x - 77)) + 6 11 
(Je * d x r> 5 (x) + Je * 8 x f e > 5 (x - r/)) sinh(r/) 



dr/ . (33) 



cosh(r?) + cos(Je * f e ' s (x) + J € * / e ' 5 (x - 77)) + 5 

For these regularized problems we show the existence of classical solutions f e ' S € C 1 ([0, T e ], H 3 (K)). 
This fact follows from the proof of the local well-posedness result in Sobolev spaces (see Section 
3.1). Now we pass to the limit in 5, showing the existence of f € E C 1 ([0, T e ], H 3 (M.)) and these 
functions are analytic. Using the previous 'a priori' estimates for i£b we conclude the existence 
of an uniform time existence, T, for all f € . Now we can pass to the limit in e, and we conclude 
the proof of the instant analyticity. 



3.3 Well-posedness for analytical initial data 

In this section we show that there is an unique local smooth solution even when the Rayleigh- 
Taylor condition is not satisfied but our initial data, which is a parametrized curve z(a), is 
analytic. For a similar result in the case with infinite depth see [6]. We prove this result by a 
Cauchy-Kowalevski Theorem (see [6, 20, 21]). 
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The complex extension of the equation can be written as 

ft* =§*•[*], 

where 

(d a z(7) - d a z(-f - r/))sinh(zi(7) - 2:1(7 - 77)) 



= P.V. f 



cosh(zi(7) - zi(j - 77)) - cos(z 2 (7) ~ ^2(7 - ??)) 



We define 



+ (dgzi(j) - d a zi(>y - n),d a z 2 (j) + d a z 2 ig - 77)) sinhQzi^) - 21(7 - 77)) 

cosh(zi(7) - 21(7 - 77)) + cos(z 2 (7) + 22(7 - 77)) ' 

j-rj/ n = sinh 2 (r//2) 

cosh(zi(7) - 21(7 - 77)) - cos(z 2 (7) ~ ^2(7 - ??)) ' 

^ + N(7,^) = w 7 s cosh 2 (77/2) 

cosh(zi(7) - zi(7 - 77)) + cos(z 2 (7) + 22(7 - V)) 

The finiteness of the term d~ is the arc-chord condition in our domain, i.e. if d" € L°° then 
z(a) 7^ z(/3) and \d a z\ > 0. The finiteness of the term d + means that the curve doesn't touch 
the boundaries. 

We consider curves z in the space 

X T = {z analytic curves on M r such that [z] < 00 and ||d + [2i]||L>x> < 00}, (37) 

where M r = {a + i£, |£| < r} and are defined in (35) and (36). We consider the following 
norm 

\\z\\ r = \\z(a) - (a,0)\\ H 3 (Mr) , (38) 

where H 3 (M r ) is the Hardy-Sobolev space (see [2] and references therein). We remark that the 
fact that they are a Banach scale can be easily proved (see [2]). For notational convenience we 
write 7 = a ± ir, 7' = a ± ir' and we take p = 2 and I = ir/2. We claim that, for < r' < r, 

C 

Pa ■ h*(B r ,) < ^—pW ' IU 2 (B r )- ( 39 ) 

Indeed, we apply Cauchy's integral formula with T = 7' + (r — r')e l ® to conclude the claim. 
Then we have the following result: 

Proposition 2. Consider < r' < r and the set 

Or = {z G X r such that \\z\\ r < R, \\d~ [z] \\L°°(M r ) < R i \\d + [z] ||i~(B r ) < R }i 

where d~[z] and d + [z] are defined in (35) and (36). Then F : Or — > X r i is continuous. More- 
over, for z,w G Or, the following inequalities holds: 

\\F[z]\\ H s (Mrl) < ^Mr, (40) 

Cr 

\\F[z] - F[w]\\ HHMrl) < ^—^\\z - w\\ H3{Mr) , (41) 

sup \F[z]fr) - F[z]fr - P)\ <C R \P\. (42) 

■yet 
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With this Proposition we can prove the local existence result: 

As zq G X ro , satisfies the arc-chord condition and does not reach the boundaries, there exists 
Ro such that zq £ O_r - We take r < tq and R > Rq in order to define Or and we consider the 
iterates 



Zn+l = zo+ F{z n ]ds, 
Jo 



and assume by induction that z^ G Or for k < n. Then, following the proofs in [6, 20, 21], we 
obtain a time Tck > 0. It remains to show that 

\\ d ~[ z n+l]\\L™(M r ) < R, \\d + [z n +l}\\L^(B r ) < R, 

for some times Ta,Tb > respectively. Then we choose T = min{TcK, Ta, Tb} and we finish 
the proof. Using the classical trigonometric formulas we obtain 

cosh (z 0l (7) -z 0l (7-?7))2sinh 2 F^zfa) - F 1 ^ - r,)ds)/2) 



+ 



+ 



sinh 2 (r//2) 

sinh (z 0l (7) - z 0l (7 - rj) ) sinh (jj F 1 [z}{^) - - ??)ds) 

sinh 2 (?7/2) 

cos (z 02 (7) " ^o 2 (7 " V)) 2 sin 2 ((£ F 2 [z]( 7 ) - F 2 [z]( 7 - 77)^/2) 

sinh 2 (77/2) 

sin (zo 2 (7) - ^o 2 (7 - ??)) sin (/<J F 2 [z](j) - F 2 [z]( 7 - rj)dsj 

sinh 2 (r//2) 



Take i < 1, assume that r] £ B(0, 1) and using the inequality (42) in Proposition 2 we have 

(d-[z rt+1 ])- 1 >-^-Ci J (t 2 + t). 

In the case where r\ E B c (0, 1), to ensure the decay at infinity, we use the inequality (40) in 
Proposition 2 to conclude 

no 

Thus we can take 



<T A < mini 1, J I — - - 



^Ro RJ 4max{C]j,C 2 } J ' 
and then 1 1 [z n +i] ||L°°(B r ) < R- We do in the same way for d + [2: n+ i] : Using the classical 
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trigonometric formulas and the previous inequalities we obtain 

cosh(z 0l ( 7 )-z 0l ( 7 -r ? ))2sinh 2 ((/ ( ;F 1 [z](7)-F 1 [z](7-r ? ) a ( S )/2) 

(d+f^!])- 1 = (d + N)- x + 



+ 



+ 



cosh 2 (77/2) 

sinh (zoi(7) - 2^(7 ~ v)) sinh (jj F 1 [z](7) - F 1 [z](-y - r])ds^J 

cosh 2 (77/2) 

cos (z 02 ( 7 ) + ^o 2 (7 " V)) 2 sin 2 ((/„' F 2 [z]( 7 ) + F 2 [z]( 7 - 77)^/2) 

cosh 2 (r//2) 

sin(z 02 ( 7 ) + z 02 (j - V )) sin (jj (7) + ^ 2 N(7 - 

cosh 2 (r7/2) 




thus we can consider 

< Tg < min 

and then ||d + [z n+ i] ||z/x>flB r ) < R- Taking < T < uuxl{Tck,Ta,Tb}, we conclude the proof. 

Remark 2 We remark that there is not any hypothesis on the Rayleigh- Taylor condition in 
this existence result. As a corollary, we obtain that if we start in the Rayleigh- Taylor unstable 
case with an analytic graph, we have local existence and uniqueness for a short time. The same 
result can be proven in the more general case of C 2,<5 functions. 



4 Differences between the two regimes 

In this section we study some properties of the regime with < A < 1 (finite depth) which are 
different with respect to the regime ^4 = (infinite depth). 

We show the qualitative behaviour for and ||9 a :/||i, 00 5 the turning waves and the ill- 

posedness in the Rayleigh- Taylor unstable case. A 'turning wave' is equivalent to have a blow 
up for d x f (see [6]). The ill-posedness result below is different from those in [9, 23] because we 
do not require a family of solutions having arbitrary long time existence. 



4.1 Maximum Principles 

In this section we show a maximum principle and a decay estimate for and a maximum 

principle for for a special class of initial data. In order to prove the maximum principles, 

the key point is to compare the local and the nonlocal terms that appear when we obtain the 
ODEs for the evolution of the L°° norms. 



4.1.1 Maximum principle for 

As in the previous section, we consider \ = I and p = 2. 

Theorem 1. Let f(t) G Hf(R) be the unique classical solution of (9) in the Rayleigh- Taylor 
stable case. Then f satisfies that 

\\f(t)h°o < ii/oiu- 
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Proof. Due to the smoothness of / in space and time we have that ||/(i)||L°° = f( x t) is Lipschitz. 
Then, using Rademacher Theorem, we have that f(xt) is differentiable almost everywhere and 
thus we get that 



d_ 
~dt 



||/(*)lk« =d t f(x t ) = P.V. [ d v f{x t - V )(E l (x t , V ,f)-E 2 (x t ,r ] ,f))dr t = I l + l2. (43) 



Let us introduce the following notation: 

a Sift) - f( x t - V) 



f(xt) + f(x t - rj) 



(44) 



Now, since Hi is defined as (10), using the classical and hyperbolic trigonometric formulas for 
the half- angle, we have 



-2P.V. / d v 9 



-2P.V. 



2 sinh (U) cosh (2) 

l (l + tan 2 (g))coth(2) 

A " 1 + tan 2 (0)coth 2 (2) 

f , 2 coth (2) 

7 M " 1 ; 1 + tan (9 coth (rj 2 2 ' 



0„0 



2 sinh (2) cosh (2) 
cosh (r?) - 1 + ^ 



-2P.V. 



d„9 coth ( 2 



an 2 (0) 
2) 



cos 2 (0)l + tan 2 (0) coth 2 (2) 



tan 



^) f ( 

- T -/- Y + P.V. / tan(0)e>„ - 



2 coth (2) 



1 + tan 2 

where we integrate by parts. Considering 

G(x) = 

we have 



+ (tan(0) coth(r//2)) 2 



1 + X' 



+ arctan(x), 



tan 



fix 



l + t an 2 (^) 



<9„ 



G 



( tan(0) 
\tanh(2J 



P.V. 



tan(0) 



sinh 2 (2)l + (tan(0)coth(2)): 



Then we get 



h 



-2f(xt) - P.V. / 

JR 



tan(0) 



-(h). 



smh 2 (77/2) 1 + (tan(0) coth(7?/2)) 2 ^ 

Anagously for the term, we use the classical and hyperbolic trigonometric formulas. In this 
case we have to write all in terms of the tan(0). This is possible because Xt is a maximum point. 
Since H2 is defined as (11) and using the same function G evaluated in tan(0) tanh(ry/2), we get 
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the following expression for I2 

f n - 2sinh(2) C osh(2) f _ 2 sinh (3) cosh (2) 

h = -2P.V. / d v 6 ^ %L = -2P.V. / d v 6 — 

7r rn.hfnl I 1 ~ tan2 ( e ) A cosh 77 - 1 + — 2 

COSH (7/ J + 1+tan 2 (e) l+tan»(0) 

(l + cot 2 (fl))coth(2) _ r C oth(f) 

" 1 + cot 2 (0) coth 2 (2) ' • 7 K s in 2 (fl) 1 + cot 2 (0) coth 2 (2) 

$,0 tanh (2) _ r tanh (§) 

cos 2 (0) 1 + tan 2 (9) tanh 2 (|) " " 2RV " 7 R " ^ } 1 + tan 2 (I) tanh 2 (2) 

= ~ 2fiXt) + L 1 + (cot(.-)coth(,/2)) 2 (46) 

Therefore, by (45) and (46), we have in (43) 



cot(fl) l_ 

cosh 2 (r//2) tanh 2 (7?/2) + cot 2 (0)' 



d t f{x t ) = -4/(x t ) + / l2 7 /oW 727 ^ , , 9/ ^ r/ 



P.V. / tan 2 i6) 1 ^dr,. (47) 

' cosh 2 (?7/2) tanh 2 (r?/2) +tan 2 (0) 



cot ((9) = tan ( ^ - 6>) = tan ( ^ - + 



Due to the definition of 

ill: = tan f , 

.2 y V2 

Moreover, using that 

00 

arctan(tan(/(x t )) tanh(r//2)) = 2f(xt), 

—00 

we can write 

J[ t} y R cosh 2 (2)tan 2 (f-/(x i ))+tanh 2 (2) a7? ' 



and we use the equality 



taU, ('I -/(„) + »)- *»» -/(*.»+*»(») 



.2 JV y / l-tan(f -/(x t ))tan(0)" 
By notational convenience we use the notation a = ^ — f(xt). We define 

„. , tan(cr) tan(0) 



tan 2 (a) + tanh 2 (2) tan 2 (6>) + tanh 2 (2) 

(tan(cr) + tan(0))(l - tan(a) tan(0)) 



(tau(o-) + tan(0)) 2 + (1 - tan(cr) tan(0)) 2 tanh 2 (2) ' 

and, using (47) we have 



dtf(xt) = - I ^2 U(x, V ,t)dr,. (48) 
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So we need to prove that II > (respectively < 0) if = maxa./(a;) (respectively 

min x f(x)). We have 

_ tan(o-)tan(0)(tan(0) + 2tan(cr)) 

~ [tan 2 (cr) + tanh 2 (2)][(tan(o-) + tan(0)) 2 + (1 - tan(a) tan(#)) 2 tanh 2 (|)] 

tan 2 (a) tan(0)(tan(<x) tan((9) - 2) tanh 2 (2) 
+ [tan 2 (a) + tanh 2 (| )][(tan(cx) + tan(6>)) 2 + (1 - tan(a) tan(0)) 2 tanh 2 (|)] 

tan(cr)tan(6»)(tan(o-) + 2tan(0)) 
[tan 2 (0) + tanh 2 (2)] [(tan(cr) + tan(0)) 2 + (1 - tan(<r) tan(0)) 2 tanh 2 (2)] 

tan(cr) tan 2 (0)(tan(cr) tan(0) - 2) tanh 2 (2) 
[tan 2 (0) + tanh 2 (2)] [(tan(a) + tan(0)) 2 + (1 - tan(cr) tan(6<)) 2 tanh 2 (|)] 

(tan(cr) + tan(0)) tan(cx) tan(0) 
+ (tan(cr) + tan(0)) 2 + (1 - tan(a) tan(0)) 2 tanh 2 (2) ' ^ ' 



+ 



Rearranging, we get 

tan(cr) tan 2 (0)[l + tan 2 (cr) tanh 2 (2)] 



n 



[tan 2 (a) + tanh 2 (2)][(tan(<r) + tan(0)) 2 + (1 - tan(a) tan(0)) 2 tanh 2 (§)] 

2tan 2 (cj) tan(0)[l -tanh 2 (§)] 



+ 



+ 



+ 



[tan 2 (a) + tanh 2 (2)][(tan(cr) + tan(0)) 2 + (1 - tan(cr) tan(0)) 2 tanh 2 (2)] 

tan 2 (cr)tan(0)[l + tan 2 (0) tanh 2 (§)] 
[tan 2 (0) + tanh 2 (2)] [(tan(cr) + tan(0)) 2 + (1 - tan(a) tan(0)) 2 tanh 2 (2)] 

2tan(a)tan 2 (0)[l - tanh 2 (§)] 
[tan 2 (0) + tanh 2 (2)] [(tan(<r) + tan(0)) 2 + (1 - tan(<r) tan(0)) 2 tanh 2 (2)] 

(tan(cr) + tan(0)) tan(cr) tan(0) 



+ 



(tan(cx) + tan(0)) 2 + (1 - tan(a) tan(0)) 2 tanh 2 (2) ' 



Now, we use that tanh 2 (2) < 1. 

If ||/||z,°o = m&x x f(x), the definitions of and a give us that tan(0), tan(cr) > and 
obtaining that II > 0. This concludes the proof for this case. For the case where the L°° norm 
is achieved in the minimum the proof is analogous. Indeed, we have that in this case II < 
because tan(cr), tan(0) < 0. □ 

Remark 3 The main difference of this result with respect to the one in [10] is that we have 
positive and negative terms in (47). Thus we have to balance them to obtain our result. We 
note that the local terms dissapear with infinite depth. 

4.1.2 A decay estimate for ||/(i)||L°° 

Theorem 2. Let /q > 0, /q € L 1 n H^(M) be the initial datum and assume p 2 — p 1 > 0. Then 
the solution f(x,t) of equation (9) satisfies the inequality 

j t \\f(t)h~ <M\\fo\yA\foh^,p\p\i)e~ J w&k. 

Proof. We conserve the notation and the hypothesis of the proof of Theorem 1, i.e. p = 2, I = 
f[x t ) = ||/(t)||x«, a = f - f{x t ) and 9 = tl&k£*t=nl m We have the equation (48) with II 
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defined in (49). Due to analysis in the proof of Theorem 1, we have the following bound 

tan 2 (V)tan(6))[l + tan 2 (0) tanh 2 (|)] 
~ [tan 2 ((9) + tanh 2 (| )][(tan(<r) + tan(#)) 2 + (1 - tan(cr) tan(#)) 2 tanh 2 (|)] 

tan(0) 1 
" (tan 2 (||/o|| Loo ) + 1)2 + tan 2 (||/o||i~) 1 + tan 2 (||/ || L -) " 

thus 

-d>f(x+)> f - tan ( g ) dri ( 50 ) 

tJ{ i; -Acosh 2 (r / /2)((tan 2 (||/o|| 1 o ) + l)2 + tan 2 (||/o|| L oo))(l + tan 2 (||/ ||L-)) 1 ' 

Fix the interval [— r, r]. We consider the sets 

Wi = j^e [-r,r],0>^f^W0 

and 

/Ml 



U 2 = yj :r] e [-r,r},6 < 4 r ■ 
The equation (13) gives us a control for the measure of these sets. Indeed, 
b||ii= / f{xt-n)df}> I f(x t - V )d V >^-\U 2 \. 



Therefore 

\U x \=<lr-\U 2 \>l{r ll/o|lil 



thus, coming to (50), 

- dt f(xA > [ - ^ dn 

tJ{ t) ~ y Wl cosh 2 (r / /2)((tan 2 (||/o|| L oo) + l)2 + tan2(||/ || ioo ))(l + tan2(||/ || Loo )) V 

> tan(/(x t )/4) 

" cosh 2 (r/2) ((tan 2 (||/ || LOO ) + 1)2 + tan 2 (||/ || L c O ))(l + tan 2 (||/ || LO o)) 1 11 

> rf(x t )-\\fo\\ L i 

" 2cosh 2 (r/2)((tan 2 (||/ || i0 o) + 1)2 + tan 2 (||/ |U~))(l + tan 2 ( 1 1 / 1 ) ) ' 

Now, we take 

r = 2 ]L/bjUi 

/(**) 

Then we get 

- dtfixt) > 



ll/ol 


L 1 


cosh 2 mf) 


|2((tan 2 (||/ || LO o) + l)2 


+ tan 2 (||/ | 


L ^))(l+tan 2 (||/ | 


Loo)) 



> 



2\\f \\ L l 

e s(xt) H/oIIli 



2((tan 2 (||/ || L ^) + 1)2 + tan 2 (||/ |U-))(l + tan 2 (||/ ||^)) ' 
and we conclude the proof. □ 

Remark 4 We observe that in the whole plane case (see [10]) the decay rate is given by 

dtf{xt) < -c(||/o||LiJ/o||L~y,p 2 )(/(z t )) 2 , 

so in the case without boundaries the decay is faster. 

Remark 5 As a corollary we conclude that there are no one-signed, integrable, steady state 
solutions. 
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4.1.3 Maximum principle for ||d:r/(i)||.L°° 

In this section we show the maximum principle for ||<9a;/(i)||.L°° for a special class of initial data: 

Theorem 3. Let fo £ Hf(M.) be a smooth initial datum such that conditions (16), (17) and (18) 
hold. Then the solution f(x, t) of equation (9) satisfies 

\\d x f(t)\\ L °° < \\d x fo\\ L °°- 

Moreover, if (x(l),y(l)) is the solution of the system (20) and assuming that ||/o||l°° < x(l) and 
\\dxfo\\L°° < y(l) } w z have that 

\\d x f{t)\\ L « < 1. 

Proof. We take p = 2 and it = 21 without loss of generality. Using the same method as in 
Theorem 1 and the smoothness of / we have that the evolution of d x f(xt) = \\d x f(t)\\L°° is 
given by 

-T7\\dxf\\L°° = d t d x f(x t ). 
at 

We suppose that d x f(xt) = m&xd x f(x,t) in order to clarify the exposition, but the proof is 
analogous in the case where the norm is achieved by the minimum. Since (9) is equivalent to 
(12), so we have to take a derivative in space in this equivalent formulation. The boundaries in 
the principal value integrals contributes with —8d x f(xt). Thus we get 



d t d x f(x t ) = -8d x f(x t )+I 1 + I 2 , 



with 



and 



t ari { iM^m) 



2RV - ' « I arCtan I tanh(^) I I *>> 
h = 2RV. J dl (arctan (tan tanh f^l 



drj. 



We define 



tanh (SfH) ' 



\ , 2 = tan (ZM±M) tanh (i 



tan iimtdM. 

\ 2 ; ( f{xt) + f{r])\ ^ , fxt-v 

I': 



We compute 



OP 2 . / >, tanh 2 ((x i - r?)/2) Q x (x t , 77, t) 

2o x arctan (/Ji x - 



with 
and 

with 



cosh 2 ({x t - n)/2) cos 2 (6») (tanh 2 ((x t - 77) /2) + tan 2 (<9)) 2 ' 



Qi = 9 x f(x t )f4 + (1 - {d x f{x t )f)iii - d x f(x t ), 

2<9 2 arctan (no) = n = Q2{xt,r], ) ^_ — 

^' cosh 2 ((x 4 -r ? )/2)cos 2 (^)(l+tanh 2 ((x 4 -r/)/2)tan 2 (0)) 2 



Q 2 = -d x f(x t )i4 + ((^/(x t )) 2 - l)/u 2 + d x f(x t ). 

Thus the sign of the integral terms are given by the sign of Q\ and Q 2 - Qi are polynomials in 
the variables /jj, respectively. 
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The roots of Q\ are d x f(xt) and —l/d x f(xt), so if we have 

1 



/ii| < min < \\d x f(t)\\ L oo, 



\dxf(t)\\ 



then we can ensure that the integral involving the increments of / is negative. However we have 
that for 7] = xt the following equality holds 



lim 

Tj-tXt 



tan 



f(xt)-f(y) 

2 



tanh (S=2) 



\dxf(t)\\L°°, 



and so we need that min{||9 x /(t)||ioo, l/\\d x f(t)\\L°o} = \\d x f(t)\\L°o. Thus we impose condition 
(16). 

Moreover, if 

I 1^1 +V, I I ^ tan (ll/o||L°°) . ,,o , n 

\x - rj\ > 1 then \/j,i\ < — — — m — < \\d x fo\\ L °o, 

tanh I | 1 

under the hypothesis (17). Then 

if = 2P.V. / d 2 | arctan j \ , ,.;_„, / | ) </// < U. 

JB°(x tl l) 

We have to bound the following integral 



tanh (2^2) 



ds/fo) (Mi) 2 + (1 ~ {d x f{x t )) 2 )in - d x f{x t ) 
B(x t ,i) sinh 2 ((a;* - r/)/2) cos 2 (6>) (l + (/^i) 2 ) 

Using the definition of 9 in (44) and the fact that \xt — 7]\ < 1, we have 

To obtain this we split as follows 

tan(0) ^ /(xt) _ tan(0) 



tanh((x t — f])/2) tanh((xj — rj)/2) 

9 ( 1 2 \ 

^ tanh((x 4 - r})/2) ^ ^ Vtanh((x 4 - rj)/2) x t - r/J Xt • 

Taylor theorem and the fact that the function (xt — rj)/ tanh(xt — rj) < 0.5/ tanh(0.5) in this 
region gives us the desired bound. In the same way, 

,an2(<,) m(„)f 



tanh 2 ((x t - rj)/2) 



(xt - n) 2 I f 9 -^ £ T A 
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Thus, using the cancellation when \i\ = d x f(xt), we obtain 



[d x f{x t ) + 5(d x f(x t )f) 1 + d x f(x t ) 8 x f(x t ) + 



tan 



( a x f(x t ) 



Til. 

1 1 



tanh(|) 



< 



48tanh(l/2)cos 2 (||/(t)|| Lo 
{d x f(x t ) + 5(d x f(x t )) 3 ) I 




B(0,1) sirm V2 

l + d x f(x t )(d x nx t ) + ^ry 



< 




We have 



h < 



Easily we get 



6tanh(l/2) cos 2 (||/(i)|| L oo) 

{(d x f{x t )) 2 - 1) tanh((x t - rj)/2) tang) + d x f(x t ) 
cosh 2 ({x t - rj)/2) cos 2 ((9)(l + tanh 2 ((x 4 - 77) /2) tan 2 (fl)) 2 



COS 2 (||/(i)|| L ~) 

It remains to show that 

if + I 2 - 8d x f(x t ) < 0. 

We need to use the local term —8d x f(xt) in order to control the remainder terms. Using the 
maximum principle (Theorem 1) we obtain 



\dxfoh^ + 5||dz/o|||o 
6tanh(l/2) 



tan 



1 + ||<9a;/o||L 



\d x fo\ 



+ 



\\Mok 



tanh(l/2) 

+ 4tan(||/ || L o O ) + 411^/011^(1 - 2cos 2 (||/ ||l-)) <0, 

which is the condition (18). 

What we have to show is that if initially the previous conditions holds there is local in 
time decay for ||<9x/(i)||L°°, but maybe these conditions are not satisfied for all time. Indeed, if 
\\9 x f(t) decays faster enough then condition (18) is not global. But if there exists ii > 
such that ||3 x /||ioo starts to grow then there is a time t* > t\ so that the condition (18) again 
holds so (19) is achieved. The same is valid for the condition (17). 

The last part in the Theorem is obvious using the maximum principle for ||/(t)||z,°° (see 
Theorem 1). This concludes the proof. □ 

Remark 6 We observe that in the case .A = the condition depends only on ||c?a;/o • In the 
case < A < 1 this appears to be impossible because of two facts: First, the term with 6 in (44) 



gives us a condition on ||/o||l°° > also we notice that the condition on tan 
gives us implicitly a condition on ||/(t)||z,°°- Indeed, we have 



2 



/ tan h(M) 



tan 



lim • 



2 



tanh 



. 2 ) 



|tan(/(z t /2))| < \\d x j{t)\\i 



The second fact is that the term \i\ can be bounded below by the incremential quotients, but if 
we want to bound it above we have to use ||<9a;/(t)||Loo and 

Remark 7 The region without decay but with an uniform bound (see Figure 2) appears due 
to the boundaries. This region does not appear in the case with infinite depth A = 0. We notice 
that if we now take the I — > oo limit we recover the well-known result (contained in [10]) for the 
whole plane case. Indeed, if / — > oo, the conditions (17) and (18) are automatically achieved 
and we only have (16) as in [10]. 
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4.2 Turning waves 

Now, we can prove the existence of turning waves in the stable Rayleigh- Taylor regime. 

Theorem 4. Take p 2 — p 1 > 0. Then there exist analytic initial data zq = z(a,0), that can 
be parametrized as a graph, such that the solution of (14), at finite time t = t\, is no longer a 
graph. 

Proof. We take p = 2 and I = tt/2. First, we show that there exists curves z(a) = (zi(a), z 2 (a)) 
such that: 

CI. Zi are analytic, odd functions. 

C2. d a zi{a) >0,Va ^ 0, d a zi(0) = 0, and d a z 2 {0) > 0. 

C3. a «i(O) = dadtz^O) < 0. 

By integration by parts in expression (14) and using the definition of Z{ we obtain 

1 



d a vi(0) = 2dz 2 (0) / d a zi(rj) sinh(^i(r?)) sin(z 2 (r?)) 
Jo 

Now, we define piecewise smooth and odd functions 

z\{a) = a — a exp 
and, fixed < b < a positive constants, 



+ 



(cosh (21(77)) - cos(z 2 (r7))) 2 
1 



(cosh(zi(ry)) + cos(z 2 (7?))) J 



dq. (51) 



z 2 la 



sin(aa) 



1 

a 

1 - / 

— sin 7T 

a V 



-a 



if < a < 



a — (ir/a) 
(Tr/a) - (tt/6) 



7T 



a 



-a + 



a 




7T + 



7T 



11 — < a < — , 
a 
. r vr 7T 

lf 6^ a< 2' , 
1 f-< a<7r( l-- 

if tt(1 - -) < a. 



(52) 



We observe that C2 is achieved for this (zi,z 2 )- Moreover, these curves satisfy the arc-chord 
condition in the whole domain. Using the definition of z 2 we have that 

d aVl (0) < 2dz 2 (0)(I a + h), 

where I a ,It are the integrals (51) on the intervals (0, vr/a) and (vr/6, tt), respectively. Easily, 
we show < and this is independent of the choice of a. The integral I a is well defined and 
positive, but goes to zero as a grows. Therefore, by approximating, there exists curves (z\,z 2 ) 
that satisfies the conditions C1-C3. 

Now, we consider (z\,z 2 ) as the analytic initial datum for the equation (14). By a Cauchy- 
Kowalevski Theorem, there exists a curve, w(a,t), solution of (14) for any t € [— T, T] (see 
Section 3.3). Due to C3, for — T < t < w(a,t) can be parametrized as a graph, at t = 
w(a, 0) = z{a) has a vertical tangent, and for < t < T the curve is no longer a graph and the 
Rayleigh- Taylor condition is not satisfied in a neighbourhood of a = 0. □ 
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Remark 8 This theorem implies that there exist initial data /q, parametrized as graphs, such 
that the solution of (9) develops a blow up for \\d x f(t) at finite time t\. 

Theorem 5. We consider a strip with width equal to I, a fixed constant. Then there exists 
initial data zo(a) = (zi(a), 22(a)) such that the solution of (14) achieve the (Rayleigh- Taylor) 
unstable case and, if you consider the same initial datum when the depth is infinite, the same 
curves are depleted. 

Proof. We can assume that p 2 — p 1 = 4-7T and take I = tt/2. It is enough to show that there exist 
curves zo(a) = (zi(a), 22(a)) such that the following holds: 

1. d a vi(0) = d a dtzi(0) > in the deep water regime, 

2. d a vi(0) = d a dtzi(0) < when the strip is considered. 

We construct a piecewise smooth curve such that both conditions holds. We take z\ defined as 
follows 



with k = le — 4. The idea is to take k « 1 such that e a k ~ 1, for — 7r < a < tt. Moreover, 
we take 22 as in (52) with a = b = 3, i.e., 



z\(a) = a — 




sin(a) 



z 2 {a) 



= < 




Notice that, in the deep water regime, the expression (51) takes the form 




Substituting the choice of z, we need to compute 



2 



(53) 



where 




4 f 3 (1 — cos(??)e r ? k + 2kr)e ??2fc sin(ry))(r/ — e r ^ k sin(r?)) sin(3r?) 
3 Jo (r? - e-^ k sm(r])) 2 + (sin(3r ? )/3) 2 ) 2 



drj 



I2 



= 4 



2 (1 — cos(n)e r,2k + 2krje r ^ k sm(rj))(rj — e r ^ k sin(T/))(— 77 + 7r/3) 
(77 - e-^ fc sin(r/)) 2 + (-r? + vr/3) 2 ) 2 



dr7, 



and 




2 



(54) 



25 



where 

X 4 = / (1 - cos(r])e~ v2k + 2krje~ v2k sin(r/)) sinh(?j - e~^ k sm(r])) sin(sin(3r/)/3) 
JO 

.( 5 1 

\(cosh(?y — e ^ fc sin(?7)) — cos(sin(3r/)/3)) 2 



1 

(cosh(?7 — e~ v i 2k sin(r/)) + cos(sin(3r/)/3)) 2 




Is= (1 - cos(r])e~ v2k + 2kr]e~ v2k sin (77)) sinh(r/ - e~' ?2 ' c sin(r/)) sin(-r? + 7r/3) 

1 

(cosh(r/ — e _r ? 2fe sin(r/)) — cos(— r/ + vr/3)) 2 

(cosh(r/ — e~^ 2fc sin(ry)) + cos(— rj + 7r/3)) 2 ^ ^' 

and 

fir 

X 6 = / (1 - cos(r?)e~' ?2fc + 2kr)e~ r]2k sinfa)) siring - e^ 2 * sin(r/)) sin (77 - 27r/3) 
■'I 

1 



(cosh(?7 — e v k sin(ij)) — cos(r/ — 2-7r/3)) 2 

+ 1 \ d 

(cosh(7y — e~ r > 2k sin(r/)) + cos(r/ — 27r/3)) 2 / 

In order to obtain the sign of (53) and (54), we compute the integrals Z$, i = 2, 3, 5, 6 using 
the trapezoidal rule with a fine enough mesh, dx = 10~ 9 (see Figure 3). The integrals Xj, i = 1, 4 
are approximated by 

, 4 /"f (1 — cos(7/)e _r?2fc + 2kr\e~ r]2k sin(ry))(r; — e _T?2fc sin(ry)) sin(3r/) 
1 = 3 io.! (77 - e-^sinfa)) 2 + (sin(3r/)/3) 2 ) 2 



and 



Z4 = / (1 - cos{ri)e' r,2k + 2A;r/e"' ?2fc sin(??))sinh(?? - e"' ?2fc sin(r/)) sin(sin(3r?)/3) 

JO.l 



(cosh(?7 — e r ' 2fe sin(?7)) — cos(sin(3r/)/3)) 2 

(cosh(r/ — e _T ? 2fc sin(r/)) + cos(sin(3^)/3)) 2 ) ^ 

We need to ensure the sign, thus we have to bound analytically the error. The truncation of the 
integral domains in XL i = 1,4 gives us an error E PV < 0.72 • 1(T 3 . The other errors are the 
one coming from the floating point representation Ep < 8.5 • 1CP 3 , and the one coming from the 
method used in the numerical quadrature Ej < 1.1 • 10~ 3 . We conclude that, if d a vi(0) denotes 
the numerical approximation of d a v\(0) defined in (54), we have 

a "i(O) < d a v\(0) + \E PV \ + \E F \ + \Ep\ < 0, 
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Figure 3: Approximating d a vi(0)/2 with different spatial step dx. 



and, in the case where d a v\(Q) is defined in (53) we get 

< d a vx(Q) - \E PV \ - \E F \ - \Ej\ < 9 a «i(0). 
Finally, we approximate this zq by analytic functions. This concludes the proof. 



□ 



This theorem shows that the problem with finite depth appears to be, in this precise sense, 
more singular than the case A = 0. 

4.3 Ill-posedness in the Rayleigh- Taylor unstable case 

Now we show that if we consider p 2 < p 1 the problem is ill-posed in Sobolev spaces (singularity 
appears for arbitrarly small initial energies and times). We remark that, if the initial data is 
analytic, there is local existence (see Section 3.3). The idea is to use the instant analyticity 
forward in time to conclude: 



Theorem 6. There exists a solution f of (9) with p 2 < p 1 such that \\fo\\ii s 
\h s (r) = 00 j f or an V s > 4, e > and small enough 5 > 0. 



< e and 



Proof. We prove the case s = 4, being analogous the rest of the cases. Take g (x) G H 3 (R) but 
go £ -ff 4 (M) and p = 2. We consider a fixed constant R > 4 and < A < 1. Now we denote 
f x (x,t) the solution to the problem (9) with initial datum f x (x,0) = Xgo(x). We know that / A 
exists for a positive time T(A, go) an d that it is analytic in a complex strip which grows with 
constant k(X,go) (see Sections 3.1, 3.2 and equation (28)). We can take an uniform k* with 
respect to A. Indeed, using the definition of (24), 



IT 



1 



7rm(0) 7r 1 ,. 

4 4 1 + \ 2 \\d x g \\ 2 Loo{R) ~ 4 1 + II^q ||| 00( 



k*(g ). 



Then the condition 



4F(<7o)-27rm(0) =4k*{g ) 



2tt 



1 + \ 2 \\d x g \\ 2 



< 
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is satisfied. We have 

l + \ 2 \\d x g \\ 2 Lo 

<Sb[/ (°)] = a ll3o||fr3( R ) + ||rf [^go]\\L™(R) + \\d + [\go]\\ L °°(R) + \\ D [^9o]\\L~>(R) + 



7T 



with 



... sinh 2 (r?/3) 2 
\d [Ago\\\L°°(UL) = = ^< 



\d + [Xg 



'Oj||L° 



2sinh i (r//2) 
cosh 2 (77/3) 



sin 2 (f (go^-gp^-??))) ^ 9' 
i sinh 2 (»)/2) 1 



2cosh 2 (r ? /2)(l- Sin2 (^^ )) ) 



2cosh 2 (?? /2) (l-^^f^ 
and ^ 

P[Aso]||l»(K) = cosh (i?) -2cosh(2A|| 5 o||L-(M)) ~ cosh(4) -2cosh(7r)' 
Therefore, we conclude the following uniform bound: 

£b[/ A (0)] < c(g ), 

with c(go) some constant depending on go that changes from line to line. We also take 

0<S*(g ) = ) < min T(X,g ). 

cexp(c(g )) Ae[o,i] 

Now, we consider < 5 < 6* (go)- We remark that all f x (x,t) exists up to time 5*(go) and, by 
means of the instant analyticity result in Section 3.3, we have 

£M[f X (t)] < c(g ), V0<t<6*(g ). (55) 

Now define f X ' 5 (x, t) = f(x, -t + 5). We have 



II/ a ' 5 (^)IIh 4 (r) = <MMI# 4 ( 



oo. 



Recall that f x is analytic in the common complex strip growing with constant k* (go) for all 
< A < 1. Then, applying Cauchy's integral formula with the curve T = x + k*(go)6e ie and 
that Hardy spaces on growing strips are a Banach scale, we get 

11^/^(0)11^) = \\d 4 J X (S)\\ LHR) < ^^—plf^^y 

Using the uniform energy bound (55) we have 

II^/ A Hl 2 (B|) < c(g )X\\d^.go\\ L 2 m . 

thus 

R 4 / A ' 5 (o)|| L2{K) < 4f±\. 

Now, given e > take < A = min jl, ^Ly j to conclude \\d*f x > s (Q) || L 2 (]R) < e. □ 
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